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•  The  difficulties  for  SOC  evaluation  of  the  LFP  cell  were  analysed. 
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•  The  modelling  of  the  voltage  hysteresis  was  introduced. 

•  Treatment  of  measurements  errors  was  done  using  a  Kalman  filter. 

•  The  algorithm  effectiveness  has  been  evaluated  through  measurements  taken  from  lab  tests. 
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This  paper  describes  a  state  of  charge  (SOC)  evaluation  algorithm  for  high  power  lithium  iron  phosphate 
cells  characterized  by  voltage  hysteresis.  The  algorithm  is  based  on  evaluating  the  parameters  of  an 
equivalent  electric  circuit  model  of  the  cell  and  then  using  a  hybrid  technique  with  adequate  treatment  of 
errors,  through  an  additional  extended  Kalman  filter  (EKF).  The  model  algorithm  has  been  validated  in 
terms  of  effectiveness  and  robustness  by  several  experimental  tests. 
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1.  Introduction 

The  lithium  iron  phosphate  (LFP)  has  emerged  as  one  of  the 
favoured  cathode  materials  for  lithium  ion  batteries,  especially  for 
use  as  an  energy  storage  device  (ESS)  in  hybrid  electric  vehicles 
(HEV)  and  electric  vehicles  (EV),  thanks  to  its  high  intrinsic  safety, 
capacity  for  fast  charging  and  long  cycle  life  [1  .  Recent  research 
and  development  in  this  technology,  especially  the  introduction  of 
the  nanoscale  LFP  cells  has  led  to  phenomenal  increases  in  its 
specific  energy.  Their  voltage  balance  of  3.3  V  per  cell  is  also  ideal 
for  new  applications  such  as  lithium  ion  starting-lighting-ignition 


List  of  acronyms:  EKF,  Extended  Kalman  filter;  ESS,  energy  storage  system;  EV, 
electric  vehicle;  HEV,  hybrid  electric  vehicle;  LFP,  lithium  iron  phosphate;  MST, 
Multiple  Step  Test;  OCV,  open  circuit  voltage;  SLI,  starting  lighting  ignition;  SOC, 
state  of  charge. 
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(SLI)  batteries  [2  .  Graphite  is  used  a  as  negative  electrode.  These 
cells  have  entered  series  production  and  are  being  marketed  as 
custom  off-the-shelf  solutions  for  the  emerging  HEV  and  EV  market 
by  a  number  of  companies  today.  The  commercially  available  high- 
power  nanoscale  LFP  cells  can  deliver  and  accept  30-35  times  their 
normal  current  rates,  and  match  the  performance  of  the  latest 
super-capacitors  over  short  time  periods. 

The  most  important  metric  for  an  electrochemical  ESS  such  as  a 
rechargeable  lithium  battery  is  the  accurate  runtime  evaluation  of 
its  state  of  charge  (SOC),  which  is  defined  as  the  percentage  of  the 
completely  extractable  charge  capacity  remaining  in  the  battery. 
The  SOC  indicates  the  amount  of  electrical  energy  remaining  in  the 
battery  pack.  An  accurate  runtime  estimate  of  the  SOC  is  important 
both  for  the  battery  application  designers  as  well  as  for  the  battery 
users.  Armed  with  the  insight  that  the  battery  SOC  would  be 
determined  accurately,  the  designer  is  able  to  efficiently  use 
available  battery  capacity  and  ensures  less  designed-redundancy  in 
its  capacity;  thus  providing  for  smaller  and  lighter  batteries.  With 


T.  Huria  et  al.  /  Journal  of  Power  Sources  249  (2014)  92—102 


93 


an  accurate  indication  of  the  battery  SOC,  the  user  ensures  that  the 
battery  is  not  over-charged  or  under-discharged;  and  suffers  less 
range  anxiety.  Overall,  the  battery  lasts  longer  and  provides  a  better 
performance. 

Over  the  years,  many  techniques  have  been  proposed  for  esti¬ 
mation  of  the  battery  SOC,  which  depend  upon  the  battery  chem¬ 
istry  and  the  final  application  [3,4].  The  most  reliable  test  for 
establishing  the  SOC  of  a  battery  is  to  charge  or  discharge  it 
completely,  thus  physically  reaching  100%  or  0%  SOC.  This  is  often 
adopted  for  an  EV  or  a  plug-in  HEV  that  are  charged  completely 
every  evening,  and  allows  the  onboard  SOC  estimation  algorithm  to 
gain  valuable  feedback  from  this  activity  to  recalibrate  itself. 
However,  for  a  usual  HEV  that  is  never  charged  from  the  grid, 
ampere  hour  counting  (also  called  the  book-keeping  system), 
which  uses  the  value  of  the  current  integral  as  a  direct  indicator  of 
SOC,  remains  the  most  popular  technique.  Since  the  integration  of 
current  offsets  or  uncompensated  and  systematic  effects  can  in¬ 
crease  the  uncertainty  of  measurements,  a  periodical  correction  of 
this  technique,  using  correction  point,  is  required. 

The  open  circuit  voltage  (OCV)  vs  SOC  correlation  curve  is  often 
used  to  provide  correction  points  for  the  recalibration,  but  it  cannot 
be  used  by  alone  for  batteries  that  exhibit  hysteresis  and  have  a 
very  flat  OCV-SOC  correlation  curve,  such  as  for  LFP  batteries. 
Moreover  these  batteries  have  capabilities  similar  to  the  super¬ 
capacitors:  thus  they  can  be  able  to  deliver  very  high  current  for 
very  short  time  periods  (almost  35  times  their  nominal  value).  For 
this  reason  small  changes  of  the  offset  or  the  effects  mentioned 
above  could  produce  an  integration  drift.  Hence,  it  is  necessary  to 
develop  a  robust  algorithm  able  to  accurately  estimate  the  instan¬ 
taneous  charge  stored  in  the  battery  with  the  characteristics 
described  above.  Different  techniques  have  been  adopted  in  the 
past  to  estimate  the  battery  SOC,  combining  them  with  hysteresis 
estimation  models  [5-7],  especially  using  the  Extended  Kalman 
filter  [8—11]. 

This  paper  presents  a  technique  to  estimate  the  SOC  of  an  LFP 
battery,  that  is  able  to  replicate  the  hysteresis  behaviour  of  the 
OCV-SOC  correlation  curve,  using  a  simple  equivalent  circuit 
model  presented  in  Ref.  [12  with  the  addition  of  an  extended 
Kalman  filter  to  reduce  the  measurements  uncertainty.  Section  2 
presents  the  challenges  for  SOC  evaluation  applied  to  LFP  batte¬ 
ries,  Section  3  explains  the  algorithm  model  definition  and  Section 
4  the  application  of  the  method  to  several  different  case  studies. 


2.  Challenges 

2.2.  Voltage  relaxation  time 

As  mentioned  in  the  Introduction,  the  OCV-SOC  correlation 
curve  is  often  used  to  correct  the  current  integral  errors  during 
runtime.  This  is  usually  done  when  the  vehicle  has  been  at  rest 
(with  its  battery  neither  charging  nor  discharging)  for  a  long  time 
(30  min-1  h),  and  when  its  battery  voltage  at  the  terminals  is 
assumed  to  approximate  the  value  of  the  OCV.  This  assumption  is 
valid  for  most  of  the  battery  chemistries.  An  attempt  was  made  to 
validate  this  assumption  for  the  nanoscale  LFP  cells  using  pulse 
discharge  and  charge  tests.  Under  this  test,  the  cell  was  first 
completely  charged,  rested  for  2  h  and  then  subjected  to  ten 
discharge  pulses  interspersed  by  1  h  rest  phases  till  the  cell  was 
completely  discharged.  Subsequently,  the  cell  was  charged  using 
ten  charge  pulses  interspersed  by  1  h  rest  phases  till  the  cell  was 
completely  charged.  The  cell  was  then  allowed  to  rest  for  13  h  Fig.  1 
presents  the  plot  of  the  cell  terminal  current  and  voltage  mea¬ 
surements  during  the  complete  pulse  discharge  and  charge 
experimental  test. 

The  primary  interest  in  the  experiment  was  in  validating  the 
assumption  that  the  voltage  of  the  nanoscale  LFP  cell  relaxes  to 
approximately  reach  its  OCV  after  a  long  rest  of  1  h.  However,  it  was 
observed  that  after  the  complete  pulse  discharge  and  charge  test 
(Fig.  1 ),  when  the  SOC  reached  100%,  the  voltage  did  not  relax  to  its 
OCV  even  after  13  h.  However,  it  was  found  that  the  voltage 
relaxation  time  was  considerably  less  between  the  30-70%  SOC 
window  (time  constant  of  around  800  s),  compared  to  when  the 
SOC  was  near  0%  or  100%. 

Additionally,  as  seen  from  the  characteristics  of  the  nanoscale 
LFP  cell  by  a  prominent  manufacturer  presented  in  Table  1,  the  cell 
has  a  current  capability  of  over  35  times  its  nominal  current  ca¬ 
pacity;  thus,  even  the  offsets  in  the  current  measurement  would 
integrate  (in  the  current  integral  SOC  evaluation  technique)  very 
rapidly,  to  reach  a  significant  fraction  of  the  total  cell  capacity  in  a 
small  time  period.  Thus,  for  an  accurate  runtime  estimation  of  the 
SOC  of  these  cells,  the  current  integral  SOC  evaluation  technique 
would  need  to  be  corrected  using  the  OCV-SOC  correlation  curve 
much  more  frequently  (at  most  every  one  or  2  h)  as  compared  to 
after  eight  to  10  h  of  operation  for  other  battery  chemistries.  This  is 
another  significant  challenge  to  be  overcome. 
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Fig.  1.  Current  and  voltage  plots  with  time  for  Multiple  Step  Test. 
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Table  1 

Test  LFP  cell  specification. 


Characteristic 

Value 

Cell  capacity  (minimum) 

4.4  Ah 

Cell  voltage  (nominal) 

3.3  V 

Energy  content  (nominal) 

14.5  Wh 

Discharge  power  (nominal) 

550  W 

Operating  temperature 

-30  °C  to  55  °C 

Storage  temperature 

-40  °C  to  60  °C 

2.2.  Hysteresis  and  flat  OCV-SOC  correlation 

As  inferred  from  the  pulse  discharge  and  charge  experimental 
test  (Fig.  1 ),  the  voltage  at  rest  was  plotted  with  respect  to  the  SOC 
for  different  times  of  rest,  to  estimate  the  marked  hysteresis  the 
nanoscale  LFP  cells  exhibit.  In  fact,  the  hysteresis  is  correlated  with 
the  long  voltage  relaxation  time,  with  the  level  of  hysteresis 
decreasing  with  the  rest  period,  as  visible  in  Fig.  2. 

This  shows  that  the  nature  of  the  hysteresis  is  not  steady  but 
changes  with  time.  Furthermore,  the  pulse  discharge  and  charge 
tests  also  produced  the  OCV-SOC  correlation  curve  for  the  cell  as 
shown  in  Fig.  3.  It  is  apparent  that  in  the  SOC  window  of  40%-60% 
the  curve  is  extremely  flat.  Moreover,  for  the  entire  length  of  the 
SOC  window  from  20%  to  80%,  the  voltage  gradient  is  of  the  order  of 
0.1  V.  Even  a  small  error  in  voltage  measurement,  would  lead  to 
extremely  inaccurate  results.  Thus,  implementing  the  usual  algo¬ 
rithm  of  using  an  OCV-SOC  correlation  curve  to  correct  the  integral 
current  to  estimate  the  SOC  of  the  battery  at  runtime  is  a  challenge. 

2.3.  Low  temperature  performance 

The  performance  of  LFP  cells  have  been  reported  in  literature  to 
be  adversely  affected  by  low  ambient  temperatures  [13-16].  Pre¬ 
liminary  tests  at  different  temperatures  were  designed  to  verify  the 
ability  of  the  cell  to  deliver  and  accept  current  at  a  required  value. 
As  seen  from  Fig.  4,  the  test  was  successful  at  20  °C.  At  0  °C,  the  cell 
demonstrated  the  required  discharge  capability,  but  was  unable  to 
charge  as  the  cell  voltage  breached  the  upper  limit  of  3.6  V.  Both 
at  -15  °C  and  -20  °C,  the  cell  was  unable  to  deliver  the  required 
power  in  both  discharge  and  charge. 

It  is  to  be  noted  that  during  charge  at  -15  °C  and  -20  °C  the 
current  did  not  reach  the  imposed  required  value,  because  of  the 


Fig.  3.  Flat  OCV-SOC  correlation  curve  for  nanoscale  LFP  cells. 


voltage  limitation.  The  charging  hardware  automatically  limited 
the  charge  current  to  avoid  cell  damage;  the  charging  phase  is  not 
immediately  stopped  so  that  the  shapes  show  what  happens.  The 
delay  between  reaching  of  cell  voltage  limitation  and  actual  stop¬ 
ping  of  the  current  is  around  half  a  second.  Thus,  the  difficulty 
encountered  in  the  nanoscale  LFP  cell  tested  was  its  inability  to 
deliver  or  accept  more  than  20%  of  the  power  specified  by  the 
manufacturer,  especially  at  temperatures  below  0  °C.  These  aspects 
require  the  presence  of  a  battery  package  thermostatically 
controlled,  to  be  continuously  maintained  at  the  reference  tem¬ 
perature  of  20  °C.  The  battery  model  was  developed  taking  into 
account  this  operating  condition. 

3.  Algorithm  formulation 

3.1.  The  battery  model 

A  typical  mathematical  model  for  modelling  electrochemical 
cells  was  shown  in  Fig.  5.  The  two  representations  are  substantially 
equivalent,  where  It  is  the  terminal  current,  i.e.  the  current  that 
enters  the  cell  terminals,  and  7m  is  the  main  current,  i.e.  the  current 
that  determines  the  charge  storage. 
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Fig.  2.  Hysteresis  decreasing  with  rest  between  pulses  in  the  Multiple  Step  Test  on  LFP  cell  (20  °C). 


T.  Huria  et  al.  /  Journal  of  Power  Sources  249  (2014)  92—102 


95 


Peak  power  test  20°C  Peak  power  test  0°C 


Peak  power  test  - 1 5°C  Peak  power  test  -20°C 


Fig.  4.  Results  of  the  peak  power  test  at  different  temperatures. 


The  left  part  of  the  lower  circuit  is  a  fiction  electric  network,  the 
capacitance  and  resistance  are  not  dimensionally  equivalent  to  the 
real  ones,  Usoc  is  a  dimensionless  voltage  between  0  and  1 
numerically  equal  to  the  cell’s  SOC  itself,  while  U0c  is  the  OCV  of  the 
cell  that  is  a  function  of  the  cell  SOC.  The  SOC  of  the  cell  can  be 
defined  as  follows: 

t 

SOC(t)  =  SOC(0)+-T  [  im(t)df  (1) 

Cq  j 
0 

where  Cq  is  the  capacity  of  the  battery  cell  (Cq  =  3600CnOm)-  The 
transient  behaviour  of  the  cell  is  modelled  by  the  network  of  re¬ 
sistances  (. Ro ,  R\  and  R2)  and  the  capacitors  (Ci  and  C2)  as  described 
in  the  set  of  equation  (2). 


/ 


SOC  = 

Usoc  1  h 

^sdOj  Cq 

RC !  - 

Urc\  ,  It 
R\Ci  +  Ct 

rc2  = 

UrC2  1  It 

R2c2  ^  c2 

Ut  =  Uoc  +  +  Urq  +  Urc2 

Additionally,  it  has  also  to  be  noted  that: 

•  The  parasitic  non-linear  resistor  Rp,  simulates  all  phenomena 

that  yield  coulombic  efficiency  lower  then  unity  for  the  cell, 

including  self-discharge. 

•  The  resistor  Rs d,  simulates  the  effects  of  self-discharge. 

The  effects  of  Rp  and  Rs d  can  be  normally  neglected  when  ana¬ 
lysing  short  time  tests  and  can  be  observed  only  when  considering 
very  long  time-span,  since  they  are  related  to  parasitic  effects, 
which  for  lithium  batteries  are  very  slow.  Indeed  it  is  known  that  if 
the  tests  are  made  lasting  short  times  (not  more  than  a  few  hours) 
the  coulombic  efficiency  for  this  type  off  cell  is  near  the  unity. 

Finally,  the  model  does  not  include  any  hysteresis  effect:  the 
numerical  values  of  voltage  sources,  resistances,  capacitances  are 
intended  to  be  function  of  SOC  and  temperature,  but  not  on  the 
current  direction. 

In  this  paper  SOC  estimation  is  based  on  an  equivalent  circuit 
model  with  generic  number  of  RC  blocks,  presented  at  the  bottom 
of  Fig.  5,  that  need  to  be  chosen.  Further  details  have  been 
explained  in  next  section. 

3.2.  Parameters  identification 

The  identification  of  the  parameters  Ro,  R\,  R2,  Ci,  C2  is  based  on 
the  experimental  Multiple  Step  Test  (Fig.  1),  that  simulates  the 
dynamic  behaviour  of  the  cell  at  different  SOC  at  constant  tem¬ 
perature  (20  °C).  Using  equation  (2),  electrical  parameters  are 


96 


T.  Huria  et  al.  /  Journal  of  Power  Sources  249  (2014)  92—102 


found  by  minimizing  the  error  between  actual  and  simulated 
voltage  profile.  The  error  function  chosen  is  of  the  type: 

£  =  \J  (^actual  —  ^model)  (3) 

The  Multiple  Step  Test  (MST)  was  performed  in  two  flavours: 
charge-based  and  discharge-based.  These  two  can  be  used  to 
determine  all  the  parameters  of  the  cell’s  equivalent  circuit,  and  the 
obtained  values  can  be  compared  so  as  to  verify  differences  in  cell’s 
behaviour  when  charging  and  discharging. 

As  already  mentioned,  the  mathematical  model  has  a  generic 
number  of  R-C  blocks,  and  the  best  value  for  n  had  to  be  chosen. 
The  decision  to  determine  the  optimum  number  of  R-C  blocks 
was  made  between  the  two  options  n  =  1  and  n  =  2,  which 
avoided  excessive  complexity  of  the  model  and  therefore  of  the 
procedure  to  evaluate  the  numerical  values  of  parameters.  To 
understand  the  process,  the  results  obtained  at  the  SOC  =  50%  are 
reported  in  Fig.  6,  which  shows  a  comparison  between  experi¬ 
mental  and  simulated  voltage  for  n  =  1  and  n  =  2  at  the  end  of  a 
discharge  pulse. 

The  comparison  shows  that  both  models  present  rather 
acceptable  results,  while  the  2-block  one  has  a  distinct  advantage. 
As  inferred  from  the  green-curve,  in  comparison  to  the  centre  of  the 
blue  experimental  band,  the  single  block  model  presents  a  final 
value  that  is  markedly  lower  than  the  actual  OCV.  Since  the  value  of 
Uoc  is  very  important  for  further  analysis  and  algorithms  for  SOC 
estimation,  it  is  advisable  to  force  this  number  the  most  reasonable, 
as  observable  in  the  experimental  data:  this  implies  that  the 
minimization  of  error  between  experimental  and  model  trends  is 
made  acting  just  on  R's  and  C’s,  while  Uoc  is  taken  as  being  known 
in  advance.  With  n  =  2,  on  the  contrary,  there  is  no  need  to 
constraint  Uoc ,  since  the  error-minimization  algorithm  automati¬ 
cally  gives  an  acceptable  value  for  this  parameter.  Moreover,  the 
two  R-C  block  model  gives  very  good  result  even  during  the  first 
part  of  the  considered  transient.  It  was  therefore  concluded  that 
two  R-C  blocks  are  important  to  keep  good  quality  of  the  model 
reproduction,  therefore  this  assumption  has  been  adopted.  If 
instead  the  charge-based  MST  is  considered,  the  result  when 
SOC  =  50%  is  shown  in  Fig.  7. 

The  corresponding  numerical  parameters,  following  the  same 
unaltered  identification  procedure,  present  different  values  during 
charge  and  discharge.  This  duality,  however,  is  a  potential  source  of 


Discharging  pulse  at  SOC  =  50% 


Fig.  6.  Voltage  behaviour  at  the  end  of  the  discharge  process.  Comparison  between 
experiment  and  model. 


Charging  pulse  at  SOC  =  50% 
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Fig.  7.  Voltage  behaviour  at  the  end  of  the  charge  process.  Comparison  between 
experiment  and  model. 


difficulties,  arising  when  the  cell  is  subjected  to  continuously 
varying  currents.  Indeed,  if  the  model  parameters  suddenly  switch 
from  charge  to  discharge  values  whenever  the  current  at  cell’s 
terminals  changes  its  sign,  it  is  not  difficult  to  foresee,  large  mis¬ 
matches  between  simulation  and  actual  values:  the  basic  idea,  that 
is  similar  to  what  often  found  in  literature  [6],  is  therefore  of 
adopting  hysteresis  only  for  Doc,  and  using  intermediate  values  for 
the  other  parameters.  Several  techniques  have  been  tested,  i.e.  the 
arithmetic  mean  evaluation  between  numerical  values  obtainable 
separately  during  end-of-charge  and  end-of-discharge,  or  mini¬ 
mizing  the  global  error  considering  both  end-of-charge  and  end-of- 
discharge  transients.  From  the  results  it  was  observed  that  the  two 
techniques  are  rather  equivalent.  The  actual  evaluation  of  numer¬ 
ical  parameter  of  the  model  can  therefore  be  done  in  either  way: 
the  choice  could  be  made  considering  only  practical  issues  related 
to  the  automatic  determination  of  the  numerical  values  from 
experimental  tests. 

3.3.  Hysteresis  model 

As  Fig.  2  shows,  the  LFP  cell  exhibit  hysteresis  from  0  to  100%  of 
SOC,  therefore  a  mathematical  model  of  the  hysteresis  behaviour 
of  the  cell  has  been  developed  17].  The  hysteresis  model  con¬ 
siders  the  cell  electrical  model  (Fig.  5)  as  a  dynamic  system  with 
respect  to  its  SOC  at  all  times;  moreover  it  emulates  the  hysteresis 
in  the  OCV  behaviour  during  both  the  major  loops  (0%-100%  SOC) 
and  the  minor  loop  (small  charge  discharge  cycle  during  a  longer 
overall  charge  or  discharge).  Indeed,  on-board  hybrid  vehicles  in 
general,  and  HCV  vehicle  in  particular,  the  SOC  window  during 
normal  operation  will  be  much  smaller  than  the  full  potential  SOC 
range. 

The  basic  idea  is  to  try  to  create  a  mathematical  model  that 
tends  to  move  Uoc  towards  the  upper  curve  when  the  cell  is 
charging,  and  towards  the  lower  one  when  it  discharges.  The 
speeds  of  this  transition  needs  to  be  evaluated  based  on  cell  tests. 

In  pursuing  this  objective,  after  some  attempts  the  following 
mathematical  formulation  has  been  proposed,  that  determines  the 
derivative  of  Uoc  as  a  function  of  the  direction  in  which  SOC  is 
moving,  and  on  how  far  actual  model  Uoc  is  from  the  target  curve. 
The  target  curve  is  the  upper  curve  (blue  in  web  version  in  Fig.  8) 
when  the  cell  is  charging,  the  red  one  (in  web  version)  when  it  is 
discharging. 
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The  algorithm  works  as  follows: 

•  when  the  time-derivative  of  Hoc  is  positive,  the  model  Hoc  is 
gradually  moved  towards  Hocch,  (the  blue  curve  in  web  version 
in  Fig.  8)  using  a  derivative  that  is  the  corresponding  target 
derivative  plus  a  correction  factor. 

•  when  the  time-derivative  of  Hoc  is  negative,  the  model  Hoc  is 
gradually  moved  towards  Hocdisch,  (the  red  curve  in  web  version 
in  Fig.  8)  using  a  derivative  that  is  the  corresponding  target 
derivative  plus  a  correction  factor. 

Therefore,  the  dimensionless  parameter  k  determines  the  speed 
of  transition  from  one  curve  to  the  other,  and  therefore  is  called 
transition  speed  of  the  algorithm. 

Some  tests  have  been  prepared  which  refers  to  the  somewhat 
idealized  voltage  behaviour  shown  in  Fig.  3.  The  algorithm  has  been 
tested  submitting  the  algorithm  (4)  to  a  sinusoidal  input  SOC  varied 
between  two  limits:  for  instance  the  result  between  SOC  =  0.2  and 
SOC  =  0.4,  using  a  value  of  k  =  10  is  shown  in  Fig.  8. 

This  figure  clearly  shows  that  when  SOC  increases,  the  model 
bloc  (green  curve  in  web  version)  tends  to  approach  the  upper  Hoc 
bound  (blue  curve  in  web  version),  while  when  SOC  decreases  it 
tends  to  approach  the  lower  bound  (red  curve  in  web  version). 

Several  other  tests  have  been  performed  to  test  the  robustness 
of  the  algorithm,  using  more  realistic  SOC  trends  able  to  simulate 
the  ordinary  vehicle  operation,  in  which  a  power  contribution  is 
also  coming  from  the  internal  combustion  engine,  that  vertically 
offsets  the  profile  in  such  a  way  that  the  charge  does  not  vary 
during  vehicle  operation  (Fig.  9). 


3.4.  Algorithm  SOC  estimator 

As  discussed  in  the  previous  section,  the  difficulties  of  the  most 
common  ways  to  estimate  cell’s  SOC  are  related  to  the  Ampere- 
hour  counting  technique,  or  to  the  OCV-SOC  correlation. 

Ampere-hour  counting  normally  creates  accumulated  errors 
that,  for  the  high  power  cells,  can  reach  10%  in  as  a  small  time  as 
20  min.  Therefore  this  technique  can  be  used  only  for  very  short 
times. 

OCV-SOC  correlation  brings  the  difficulty  that  for  the  central 
SOC  zone  OCV  is  nearly  horizontal,  and  that  there  is  a  marked 
hysteresis  phenomenon  that  further  complicates  the  analysis. 

The  basic  idea  that  was  issued  is  to  cycle  between  ampere- 
hour  counting  and  OCV  estimation,  to  compensate  the  errors. 
The  algorithm  basic  arrangement  can  therefore  be  as  shown  in 
Fig.  10. 

The  measures  of  the  cell’s  voltage  and  current  are  taken  from  the 
physical  system,  as  shown  in  the  upper  part  of  the  figure.  These  two 
measures  are  fed  into  the  Flybrid  SOC  evaluator,  that  then  gives  the 
output.  The  evaluator  operates  in  two  directions: 

•  In  the  upper  chain  it  makes  a  voltage-based  SOC  evaluation: 
from  the  measured  cell’s  terminal  voltage  Ht,  determines  the 
open-circuit  Hoc  (by  subtraction  of  the  voltage  drop  across  the 
electric  network  R0,  R^C-[,  R2C2)  and  then,  using  the  above  dis¬ 
cussed  SOC-OCV  correlation,  estimates  the  SOCJJ. 

•  In  the  block  containing  l/(Cs)  (“in  which  “s”  is  the  Laplace  var¬ 
iable)  another  SOC  estimation  is  obtained,  based  on  ampere- 
hour  counting. 


Fig.  8.  Effect  of  the  algorithm  in  a  sample  case. 
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The  two  SOC  evaluations  must  then  be  somehow  mixed  to 
determine  the  final  estimate  SOC_ev.  In  addition  to  the  technique 
depicted  in  Fig.  11,  adequate  treatment  of  measuring  errors  must  be 
made,  for  instance  by  taking  advantage  of  Kalman  filters.  After 
some  attempts,  the  final  implementation  of  the  basic  logic  shown 
in  Fig.  10,  which  showed  a  good  estimation  quality  as  represented 
in  Fig.  11  (where  the  upper  part  of  Fig.  10  is  omitted  for  the  sake  of 
simplicity),  was  obtained. 

In  comparison  to  Fig.  10,  the  following  differences  have  been 
introduced:  instead  of  a  simple  hysteresis  block,  a  full  Extended 
Kalman  filter  has  been  introduced,  that  performs  SOC  evaluation 
taking  also  into  account  the  measured  current  Jm,  (Fig.  5)  and 
adequately  treating  measurement  errors. 

The  Extended  Kalman  Filter  (EKF)  is  an  optimum  state  estimator 
for  nonlinear  systems,  where  a  linearization  process  is  applied  to 
the  OCV-SOC  correlation  curve  at  each  step,  in  order  to  approxi¬ 
mate  the  nonlinear  system  as  a  linear  time  varying  (LTV)  system. 
The  EKF  is  able  to  estimate  the  best  possible  values  of  the  output 
based  on  input  containing  unmeasured  noise. 

The  set  of  equation  (2)  describes  the  dynamic  linear  system 
represented  in  Fig.  5,  reproducing  the  battery  transient  behaviour 
whose  parameters  are  functions  of  the  SOC  and  changing  with  time. 
Considering  equation  (2),  it  is  possible  to  write  the  discrete  form  of 
the  battery  transient  equations: 


analysis  performed  in  the  previous  section  and  summarized  by 
equation  (4). 

Moreover,  w/<  and  v/<  are  the  process  and  the  measurement 
noises  respectively.  They  are  assumed  to  be  independent  to  each 
other  and  white  normal  distributions  with  a  zero  mean  and 
covariance  of  known  values.  Further  details  about  Kalman  filter  are 
reported  also  in  Refs.  [18,19]. 

The  other  difference  is  that  the  “mix”  block  is  chosen  as  a  switch 
that  periodically  resets  the  value  of  SOC_Ah :  normally  is 
SOC_ev  =  SOC_Ah ;  however  periodically  it  is  forced 
SOC_ev  =  SOCJcf,  and  the  integral  l/(Cs)  is  reset  to  accommodate 
for  this  new  output  value. 

In  Section  2.1  it  was  seen  that  for  high  power  cells  Ampere-hour 
counting  leads  to  errors  on  SOC  of  about  10%  after  20  min.  Therefore 
the  periodicity  of  the  mix  block  must  be  lower  than  this  value. 
Several  simulations  were  performed  on  this  time,  and  it  was  seen 
that  values  between  1  and  10  min  give  the  best  results. 

4.  Tests  and  results 

Tests  and  results  are  reported  in  reference  to  the  tests  made  on 
the  new  and  on  the  aged  LFP  cell.  During  the  tests  the  SOC  evalu¬ 
ation  is  analysed  both  for  few  hours  duration  and  long  term  sta¬ 
bility,  that  required  the  SOC  evaluation  to  be  extended  for  tens  of 
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where  Ro ,  Ri,  R2 ,  C\  and  C2  are  the  cell  parameters  as  a  function  of 
SOC  and  thus  indirectly  varying  with  time.  The  vector  of  state 
variables  x  is  constituted  by  the  two  voltages  across  R-C  blocks. 

As  usual,  matrix  A  represents  the  dynamic  evolution  of  state  x, 
matrix  B  indicates  the  quota  of  input  directly  transferred  to  states 
(algebraic  part  of  R-C  blocks)  and  matrices  C  and  D  indicate  the 
influence  of  state  and  input  into  output,  constituted  by  the  battery’s 
voltage  at  its  terminals. 

The  Extended  Kalman  filter  technique  is  applied  to  the  dynamic 
system  that  describes  the  transient  behaviour  of  the  left  part  of  the 
electrical  circuit  illustrated  in  Fig.  5.  The  EKF  uses  the  measured 
current  (7t)  and  the  OCV  (obtained  as  the  difference  between  the 
output  of  the  R-RC-RC  Network  block  and  the  voltage  measure¬ 
ments)  to  perform  the  runtime  SOC  evaluation.  The  detailed 
equations  of  this  block  are  as  follows: 


hours.  For  the  tests  on  the  standard  LFP  cell,  they  are  conducted 
simulating  the  battery  stress  on  mainly  two  reference  road  con¬ 
ventional  cycles.  After  the  identification  of  the  cell  parameters,  the 
EKF  algorithm  is  applied.  The  tests  are  repeated  taking  as  reference 
the  same  cycles,  after  several  months,  also  for  the  aged  cell. 

4.1.  Test  definition 

The  quality  of  the  SOC  estimation  applied  to  the  cell  is  tested  on 
the  repetition  of  two  main  reference  cycles.  The  two  cycles,  cor¬ 
responding  to  the  real  usage  of  the  battery  on-board  hybrid  vehi¬ 
cles,  are  shown  in  Fig.  12.  Currents  are  negative  when  the  battery  is 
delivering  power,  positive  otherwise.  Both  of  the  cycles  are  scaled 
at  cell  level,  and  a  charge  compensation  is  introduced. 

The  tests  are  executed  by  the  repletion  of  several  cycles,  per¬ 
forming  the  SOC  evaluation  through  the  algorithm. 
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Note  that  g(  )  is  the  Us oc~Uoc  correlation  function,  that  is 
able  to  take  into  account  voltage  hysteresis,  according  to  the 


4.2.  Tests  on  the  LFP  cell 

The  test  is  applied  on  the  LFP  cell  with  reference  to  the  Cycle  1. 
The  cell  is  subjected  to  the  current  profile,  then  the  SOC  estimator  is 
applied. 

The  reference  SOC  is  obtained  by  numerical  integration  of  cur¬ 
rent,  measured  with  lab  instrumentation,  and  the  parameters  of  RC 
network  are  determined  as  a  function  of  SOC,  then  the  SOC  esti¬ 
mation  is  performed  with  the  algorithm  and  finally  compared  to 
the  measured  value.  It  must  also  be  noted  that  the  proposed  al¬ 
gorithm  should  be  able  to  be  implemented  onboard  real  buses, 
using  commercial,  inexpensive  measuring  hardware.  To  show  how 
the  algorithm  is  able  to  perform  when  fed  with  measures  from 
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inexpensive  sensors,  the  voltage  and  current  data  supplied  to  it  are 
altered  with  respect  to  the  values  measured  with  lab-grade 
instrumentation  as  described  in  equation  (7).  Current  is  measured 
using  a  shunt,  class  0.5  and  full  scale  150  A;  the  voltage  is  measured 
with  a  DAQdevice  (Model  NI 9219)  with  a  declared  accuracy  of  0.3% 
of  reading. 

/  Acse/^  =  1.03/LAB(t)  +  0.2 

I  W*  =  1-02/lab (t)  4  0.04  l/J 

i.e.,  a  200  mA  offset  and  3%  deviation  is  introduced  onto  the  current, 
and  a  40  mV  offset  and  2%  error  is  added  to  the  voltage.  Moreover 


the  initial  SOC  is  set  at  70%,  although  the  real  initial  SOC  of  the  cell  is 
about  50%.  The  test  lasted  45  h.  Fig.  13(a)  shows  results  related  to 
the  first  part  of  the  test,  while  some  comments  about  long-term 
stability  are  reported  later.  The  estimated  SOC  by  the  algorithm 
(red  curve)  matches  the  SOC  obtained  by  the  numerical  integration 
of  the  experimental  measured  current  (blue  curve),  recovering  the 
initial  error  in  a  few  hours. 

4.3.  LFP  aged  cell:  results 

After  the  tests  previously  described,  the  cell  has  one  year  rest. 
After  this  period,  Cycle  2  is  applied  to  the  aged  cell.  At  first  level,  to 
evaluate  the  stability  of  the  algorithm,  no  modification  for  the 
electrical  parameters  (%  Ru  R2,  Ci,  C2)  is  applied.  At  the  same,  also 
the  EKF  Wk  and  v/<  parameters  that  influence  treatment  of  errors,  are 
kept  unmodified.  These  parameters  are  normally  tuned  off-line  for 
each  set  of  electrical  parameters  (%  Ru  R 2,  Ci,  C2)  related  to  the 
specific  cell.  Finally,  also  voltage  and  current  errors  provided  to  the 
experimental  data  fed  to  the  algorithm  present  the  same  trend  as 
described  in  equation  (7). 

The  test  lasts  about  45  h  and  is  divided  in  three  parts:  a  first 
sequence  of  48  cycles,  lasting  about  5  h,  a  second  phase  of  rest 
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Fig.  10.  Hybrid  SOC  evaluator,  basic  concept. 


Fig.  11.  Hybrid  SOC  evaluator,  actual  implementation. 
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(35  h)  and  the  second  sequence  (5  h).  Also  for  the  previous  case 
results  are  related  to  the  first  part  of  the  test,  while  some  comments 
about  long-term  stability  are  reported  later.  The  estimated  SOC  (red 
curve)  matches  the  SOC  obtained  by  the  numerical  integration  of 
the  measured  current  (blue  curve)  as  visible  in  Fig.  13(b),  recovering 
the  initial  error  in  the  same  time  used  in  the  previous  test  of 


Fig.  13(a).  This  aspect  is  representative  of  the  high  level  of  stability 
and  reliability  of  the  proposed  algorithm,  able  to  work  properly  also 
without  redefinition  of  the  parameters. 

The  second  test  is  oriented  to  evaluate  the  algorithm  through  a 
re-calibration  of  the  electrical  and  the  EKF  parameters:  particu¬ 
larly,  a  new  set  of  electrical  parameters  ( Rq ,  R\,  R2 ,  Ci,  C2)  and  the 


time  (h) 


Fig.  13.  Experimental  evaluation  of  battery  SOC  estimation  using  the  EKF-based  model  for  Cycle  1  (a),  for  Cycle  2  using  an  aged  cell  (b),  comparison  between  results  obtained  both 
for  the  aged  cell  with  original  (red)  and  renewed  (green)  set  of  electrical  and  EKF  parameters  (c).  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  article.) 


T.  Huria  et  al.  /  Journal  of  Power  Sources  249  (2014)  92—102 


101 


O 

O 

oo 


0.75 
0.7 
0.65 
0.6 
0.55 
0.5 
0.45 

0  5  10  15  20  25  30  35  40  45 


Measured  SOC 
EKF  SOC -Original 
EKF  SOC  -  Renewed 


time  (h) 


Fig.  14.  Experimental  validation  of  the  EKF-based  model  on  the  Cycle  1  (a),  validation  of  the  EKF-based  model  on  the  Cycle  2  (b). 


SOC-OCV  correlation  curve  are  estimated.  In  this  case,  providing 
the  same  voltage  and  current  errors  to  the  experimental  data,  the 
results  carried  out  by  the  algorithm  are  less  accurate  than  ones 
showed  in  Fig.  13(a).  In  order  to  have  the  same  behaviour  of  the 
last  test  SOC,  also  the  EKF  parameters  w/<  and  v/<  have  to  be 
changed.  Final  result  is  reported  in  Fig.  13(c).  This  aspect  shows 
how  the  re-calibration  of  the  electrical  parameters  should  be 
necessarily  accompanied  with  the  update  of  the  EKF  parameters, 
to  maintain  good  quality  for  the  SOC  estimation.  The  procedure, 
fully  automatable,  can  also  be  performed  during  the  standard 
vehicle  maintenance  operations,  around  once  a  year. 

4.4.  Long  term  stability 

The  long  term  stability  for  the  algorithm  is  analysed  in  reference 
to  the  three  case  studies  already  presented  above.  It  is  possible  to 
observe  that  results  already  mentioned  for  the  LFP  cell,  subjected  to 
the  Cycle  1,  show  the  stability  of  the  algorithm  for  the  total  duration 
of  the  test,  of  about  45  h.  Results  are  reported  in  Fig.  14(a). 

Similar  results  are  obtained  for  the  aged  cell  subjected  to  Cycle  2 
without  any  electrical  or  EKF  parameters  recalibration.  It  must  be 
also  noted  that  with  respect  to  the  previous  analysis,  the  SOC 
evaluation  is  stopped  during  the  rest  phase  of  the  vehicle:  also  in 
this  case  the  stability  seems  to  be  guaranteed.  Finally,  redefining 
the  electrical  and  the  EKF  parameters  for  the  aged  cell,  the  evalu¬ 
ated  SOC  matches  the  SOC  measured  through  lab  instruments  more 
accurately,  as  reported  in  Fig.  14(b). 

5.  Conclusions 

This  paper  has  mainly  put  in  evidence  the  huge  difficulties  for 
SOC  evaluation  of  the  LFP  cell’s  arising  from  its  characteristics.  First 
of  all  the  high  power,  that  enlarges  the  accumulation  effect  of 
measured  errors  on  the  ampere-hour  counting.  Furthermore,  the 
flat  OCV-SOC  curve  for  most  practical  ranges  of  SOC,  that  causes 
difficulties  in  evaluating  SOC  from  open-circuit  voltage  measures, 


combined  with  the  high  level  of  hysteresis  of  the  OCV-SOC  curve 
correlation. 

These  difficulties  have  been  dealt  at  first  modelling  the  voltage 
hysteresis:  transition  between  charge-based  and  discharge-based 
curves  has  a  higher  derivative,  than  the  stabilized  OCV  curves, 
that  eases  SOC  evaluation.  Additionally,  an  adequate  treatment  of 
measurement  errors  was  done  using  an  Enhanced  Kalman  Filter. 
The  resulting  algorithm  mixes  ampere-hour  counting,  OCV-SOC 
evaluation,  with  an  additional  correction  by  a  Kalman  filter. 

The  algorithm  effectiveness  has  been  evaluated  firstly  through 
measures  that  were  taken  from  lab  tests  with  the  highest  accuracy. 
Then  the  integral  of  these  high  accuracy  currents  was  taken  as  the 
best  estimate  of  cell’s  SOC.  Finally,  the  algorithm  has  been  fed  with 
data  that  were  altered  so  that  it  was  able  to  simulate  measures  that 
could  be  taken  from  cheap,  industrial  sensors.  Moreover  a  wrong 
initial  SOC  estimate  was  supplied  to  it. 

Evaluation  was  made  in  several  cases,  and  conditions.  In  all 
cases  the  algorithm  recovered  the  initial  errors  in  reasonable  time, 
and  remained  stable  very  near  to  the  “actual”  SOC  value.  The  al¬ 
gorithm  has  also  shown  “long  term  stability”.  Naturally,  when  the 
vehicle  is  at  rest,  e.g.  during  nights,  it  must  be  made  aware  of  this 
rest,  otherwise  it  could  continuously  crunch  just  measurement 
errors  and  offsets,  thus  leading  to  unacceptable  results.  Finally,  the 
algorithm  gave  good  results  even  after  one  year  of  calendar  life  had 
passed,  which  means  that  no  important  drift  in  the  cell’s  behaviour 
had  occurred.  However,  it  is  recommended  that,  when  the  cells  and 
batteries  are  used  onboard  production  vehicles,  yearly  recalibration 
is  made,  to  keep  the  algorithm  in  touch  with  the  actual  state  of  life 
of  cell  and  battery. 
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